Zeeman Splitting of Photonic Angular Momentum States in Gyromagnetic Cylinder 
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We show that under the presence of a static magnetic field the photon eigen- frequencies of a 
circular gyromagnetic cylinder experience a splitting that is proportional to the angular momentum 
density of light at the cylinder surface. Such a splitting of the photonic states is similar to the Zeeman 
splitting of electronic states in atoms. This leads to some unusual decoupling properties of these 
non-degenerate photonic angular momentum states, which are demonstrated through numerical 
simulations. 
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Recently, the effect of static magnetic field on photonic 
states has attracted a lot of attention due to the discovery 
of protected photonic chiral edge states in gyromagnetic 
photonic crystals [H-Qj which are analogue to the topo- 
logically protected edge states in electronic systems [J]. 
However, analytical studies on simple photonic states are 
still rare, compared to the extensive studies of electronic 
states. One of the most fundamental examples for elec- 
tronic states is the famous Zeeman effect 0. In analogy 
with the Zeeman splitting of electronic states in atom, 
here we study a similar basic effect on photonic states in a 
gyromagnetic cylinder. We also introduce some unusual 
wave decoupling phenomena which are consequences of 
such a splitting. 

Since Zeeman splitting of electronic states is associ- 
ated with the broken degeneracies of electron states with 
different angular momenta, we expect to see a relation 
between the Zeeman-like effect in our photonic system 
and the angular momentum of light. Angular momen- 
tum of light Q , which could be used for storage of quan- 
tum information Q, has drawn serious interest in recent 
years. Therefore, it is our motivation in this Letter to 
calculate the frequency splitting of photonic angular mo- 
mentum states and derive a formula which relates the 
angular momentum of light to the frequency splitting. 

We consider the splitting of photonic states in gyro- 
magnetic cylinder as shown in Fig. 1. A gyromagnetic 
material can have a permeability depending on the ap- 
plied static magnetic field. The magnetic permeability 
tensor can be written as 
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FIG. 1: (color online) Zeeman splitting of photonic states in 
the presence of static magnetic field. The top figure shows 
a gyromagnetic cylinder surrounded by air. Formula in the 
dashed box displays the relation between frequency shift and 
angular momentum density of light (j z ). The first column 
(black lines) shows the original degenerate states in the gyro- 
magnetic cylinder in the absence of magnetic field Ho — 
(i.e., with isotropic permeability /i r ). The second column 
(green lines) shows the shifts due to the change of effective 
permeability from fi r to fi' r (indicated as "index shift" ) when a 
static magnetic field of Ho = 800 Oe is considered. The third 
column (red and blue lines) shows the final Zeeman splitting 
with the adjacent signs indicating the sign of m. Only the 
states with n = 1, 2 and |m| = 4, 5 are shown. 
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For transverse electric (TE) polarization with electric 
field along the z direction, the frequencies (uj) of the pho- 
tonic states, can be calculated from the roots of the Mie 
resonance condition [8f, 
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where m is the azimuthal quantum number, J m is the 
?nth-order Bessel function, Hm^ is the mth-order Hankcl 
function of the first kind, k' = k$^J fi' r ed, fco = w/c is free- 
space wave number, and fi r = (/i 2 . — /i 2 )//i r . It should 
be noted that Eq. ([2]) does not have real root frequencies 
because of radiation loss but we will only consider the 
states with relatively low radiation loss (i.e., frequencies 
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with small imaginary parts). By comparing Eq. ([2]) with 
the Mie condition for a cylinder with isotropic permeabil- 
ity \i r [i.e., Eq. ^ with /z^ = 0], we can interpret the 
effect of the static magnetic field as two steps: (i) a shift 
associated with an index change (from permeability \x r 
to fjb' r ) and (ii) a splitting of frequencies. 

As an example for numerical demonstration, we con- 
sider yttrium- iron-garnet (YIG), which is a type of com- 
mercially available gyromagnetic materials, as the ma- 
terial of the cylinder supporting photonic states at mi- 
crowave frequencies. For an applied static magnetic field 
Hq in the z direction, we have u r = 1 +u m W(,/(^ — lu 2 ) 
and Hk = w m w/(^ — oj 2 ) [9j, where ujh = ^Hq is 
the precession frequency, 7 is the gyromagnetic ratio, 
u> m = 4njM s , AttM s is the saturation magnetization. 
Using parameters provided in a previous experimental 
study (e d = 15.26, H = 800 Oc, and 4ttM s = 1884 G) 
[lpj . we plot the frequency splitting diagram in Fig. 1 for 
a cylinder of radius a — 1 cm by finding the frequency 
roots of Eq. ^ . Here, in addition to the azimuthal quan- 
tum number m, we denotes n as the quantum number in 
the radial direction and (n, m) as a specific photon state. 
To have a clear picture, we first focus on \m\ = 4,5 in 
the lowest (n = 1) and the second lowest (n = 2) radi- 
ally quantized levels. In the absence of the static mag- 
netic field (first column of Fig. 1), the positive-m and 
negative-m states are degenerate. However, when the 
static magnetic field is present, the original degenerate 
states shift up to higher frequencies (second column of 
Fig. 1) and split into two counter rotating states (n, m) 
and (n, —m), as indicated respectively by the red lines 
and blue lines in the third column of Fig. 1. Physically, 
the effect of the static magnetic field can be understood 
as a broken time reversal symmetry (and reciprocity) so 
that photonic states are no longer degenerate. We will 
show in a later part of this Letter that such a splitting is 
proportional to the angular momentum of light. 

Using a finite element solver (COMSOL Multiphysics) , 
we examine the electric field profile of the non-degenerate 
rotational states after splitting. Fig. 2(a) and (b) show 
the wave profiles excited by a line current source ly- 
ing along the ^-direction on the surface of gyromagnetic 
cylinder at frequencies / = 8.72 GHz and / = 10.40 GHz, 
respectively. In Fig. 2(a), we see a clockwise rotating 
field pattern with four complete oscillations along the 
azimuthal directions, which shows that only the (1,-4) 
state can be excited at / = 8.72 GHz. The field pat- 
tern at / = 10.40 GHz [Fig. 2(b)] is similar to that in 
Fig. 2(a) except that it is for (1,4) state with its electric 
field rotating in anti-clockwise direction. These results 
confirm the Zceman splitting diagram in Fig. 1. 

For a complete picture of the splitting of photonic 
states, we plot in Fig. 2(c) the exact solutions [roots of 
Eq. fl2J] for the resonant frequency as a function of the 
azimuthal momentum number m. It is found that the res- 
onant frequencies show pronounced differences between 
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FIG. 2: (color online) (a) and (b) Electric field patterns ex- 
cited by a line source in the presence of external magnetic 
field (Ho = 800 Oe). (a) Only m = —4 mode is excited at its 
resonant frequency / = 8.72 GHz. (b) Same as (a), but for 
m = +4 at its resonant frequency / = 10.40 GHz. (c) Reso- 
nant frequency versus azimuthal momentum number ra. The 
solid lines with open circles show exact results obtained from 
Mie resonant conditions [Eq. Q]. The dashed lines with solid 
circles indicate approximated results obtained from the per- 
turbation theory [Eq. @]. Blue lines and red lines indicate 
the negative m modes and the positive m modes, respectively. 
Low-quality states with \m\ < 4 are not shown. 



positive and negative m states, stemmed from the effect 
of the static magnetic field. For n = 1, the most obvi- 
ous frequency splitting can be observed and the splitting 
gap width decreases gradually as m increases. We also 
observe that the splitting is relatively weak for larger n 
and, for each fixed n, the splitting is almost independent 
of m when m is large. 

To understand the frequency splitting and its relation 
to the angular momentum of light, we employ a Hamil- 
tonian eigenvalue approach for electromagnetic waves 



111 Il2| . In this approach, we construct the following 



eigenvalue problem: 
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e(r) and /2(r) are, respectively, the position dependent 
relative permittivity and permeability tensor, E(r) and 
H(r) are, respectively, the dynamic electric and magnetic 
fields. Here, we take the approximation that the disper- 
sive material properties is weak in the frequency range 
where uj is much larger than Uh ■ For TE polarization, we 
can simplify Eq. ^ as 
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and is the unit step function with a value of one 
(zero) when x > (x <0). The operator f2 is the corre- 
sponding frequency operator for materials with isotropic 
permeability ^(r) and permittivity e(r), where fi(r) = fx' r 
and e(r) = for r < a and /x(r) = e(r) = 1 for r > a. 
For small <C fJ- r , we can consider the second term in- 
volving AS1 as a first order perturbation to fin,- By doing 
so, the effect of the static magnetic field is again seen as 
the two-step process shown in Fig. 1. 

The frequency shift of each state in the first-order per- 
turbation theory is Aoj nm = -f^ 2 (ipo,n,m\AQ\tp 

0,n,m/ 

(where \ipo, n ,m) is the unperturbed state ket) and can 
be simplified as (see supplementary materials for deriva- 
tions) 
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^j-\E z {k' n a)\ 2 is the "angular momen- 



where j z \ r=a - 2u 
turn density" evaluated at the surface of the cylinder, uj' n 
is the frequency of the unperturbed state, k' n = oj' n /c, 
and U = | f{e Q e{r)\E\ 2 + n ^(r)\H\ 2 ]dA. Here, the "an- 
gular momentum density" at the boundary is defined as 
j = ^-r x [E* x (V x E)] Q. In the non-dispersive 
approximation, U can be considered as the total elec- 
tromagnetic energy integrated on the xy-plane (see sup- 
plementary materials for a discussion on the region of 
integration). Since we consider frequencies far above ujf,,, 
the third factor in Eq. ^k/(nl — n 2 .), is positive and 
the sign of frequency shift is, therefore, the same as the 
sign of m. We thus conclude that the frequency shift is 
proportional to the "angular momentum density" at the 
cylinder surface per photon. 

In comparison, the Zeeman effect on electronic states 
in atoms has a similar formula for the energy shift: 
AE = fiBgJzB/h, where J z is the total projected an- 
gular momentum in the direction of the static magnetic 
field (-B), hb is the Bohr magneton, h = h/2ir is the re- 
duced Plank's constant, and g is a dimensionless constant 
which depends on the type of atoms. Such an energy 
shift is proportional to the magnetic moments associated 
with the angular momentum of electrons. Here in Eq. 
(JSj> , we also have the proportionality between frequency 
shift and angular momentum of light. With such pro- 
portionality, wc thus call our work the Zeeman splitting 



of photonic angular momentum states in gyromagnctic 
cylinder, which is the most important result in this Let- 
ter. It should be noted that Eq. ((8]) are not limited to 
microwaves. One could easily extend this new theory to 
terahertz or optical frequencies in other systems, such as 
plasma systems where the roles of electric and magnetic 
fields are switched. 

To verify our analytical formula [Eq. (|8])]. we evalu- 
ate and compare it with the exact results in Fig. 2(c). 
Using the fact that the radiation field outside the gyro- 
magnetic cylinder should not contribute significantly to 
the frequency shift, we get a closed- form expression for 
Eq. ([8]) (see supplementary materials for derivations): 
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As shown in Fig. 2(c), we have a very good agreement 
between the frequencies evaluated from the closed-form 
solution [Eq. (j9|)] and the exact results given by Eq. ((2|) . 
Discrepancies exist only in the first few low-order reso- 
nances, which may be due to the first-order perturbation, 
the neglected radiation field, and the dispersive property 
of gyromagnctic materials. This proves the validity of 
Eq. (0). 
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FIG. 3: (color online) Proposed experimental measurement 
for the Zeeman-like splitting of photonic angular momentum 
states. The cylinder in Fig. 1 is placed between two straight 
waveguides with a refractive index 2.2 and width 5.6 mm. 
The gap between the cylinder and each waveguide is 1.8 mm. 
Arrows indicate the sources of the incoming waves, (a) and 

(b) Static magnetic field is absent. Waves launched from the 
left (a) and right (b) can also couple to the other waveguide. 

(c) and (d) Static magnetic field is present. Wave launched 
from the left (c) can couple to the other waveguide but that 
from the right (d) cannot. 

With a non-dcgcncratc angular momentum state that 
allows a rotation of electromagnetic field in only one di- 
rection, the Zeeman-like effect can lead to very unusual 
mode decoupling properties that could be useful for new 
photonic applications. For example, we consider the 
setup shown in Fig. 3. In the usual case, if we place the 
resonant cylinder between two waveguides, it can switch 
the light path from one waveguide to another waveguide 
through evanescent coupling. As shown in Fig. 3(a) and 
(b) for the cases without the static magnetic field, the 
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device is a two-way coupler because the resonator sup- 
ports both degenerate modes which couple the guided 
waves no matter which direction the source waves come 
from [see Fig. 3(a) for left-incident wave and Fig. 3(b) 
for right-incident waves]. However, by using the non- 
degenerate state (1,-4), we find that the left- incident 
wave can couple to the other waveguide [Fig. 3(c)] but 
the right-incident wave fails to interact with the cylin- 
der and pass straightly through the top waveguide [Fig. 
3(d)]. The one-way transport properties can be under- 
stood by the mode coupling [Fig. 3(c)] and mode decou- 
pling [Fig. 3(d)] between the incident guided wave and 
the single-mode state. The unusual properties shown in 
Fig. 3 can be used for experimental demonstration of the 
Zccman-likc splitting. Recently, the nonreciprocal trans- 
port properties of light in the presence of static magnetic 
field has been widely studied H3, EMU . The phe- 
nomena here may also be particularly useful for designing 
special photonic devices such as one-way coupled waveg- 
uides, diode for photonic circuit, or some add-drop filter 
systems 1_7-2(| with a much smaller region of applied 



static magnetic field (an area of one cylinder) 
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anti-clockwise rotation because this state is not allowed 
at the same frequency. This forces the electromagnetic 
energy almost localized in the first cylinder next to the 
source even when we tune the excitation frequency to 
the resonant frequency of these identical cylinders [Fig. 
4(b)]. Such a great suppression of the interaction be- 
tween resonators may be used for the design of slow-light 
devices and some applications where a high density of 
independent resonators is needed. 

In conclusion, we showed analytically that the fre- 
quency splitting of the angular momentum states in a 
circular gyromagnetic cylinder under the presence of a 
static magnetic field is proportional to the angular mo- 
mentum density of light at the surface of the gyromag- 
netic cylinder. These photonic angular momentum states 
give some unusual decoupling properties, which could be 
useful for designing novel photonic devices. Our study 
will benefit to the understanding of the angular momen- 
tum of light. 
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ported by China Scholarship Council and Southeast Uni- 
versity. We thank Anshuman Kumar, Dr. Zheng Wang, 
and Dr. Jack Ng for useful discussions. 




FIG. 4: (color online) Decoupling among photonic angular 
momentum states in an array of closely-packed identical gy- 
romagnetic cylinders. The white arrows indicate the location 
of the point source, (a) Strong coupling when static magnetic 
field is absent, (b) Weak coupling when static magnetic field 
is present. 

Another consequence of the Zeeman-like splitting is a 
localization due to the decoupling among closely spaced 
cylinders. As shown in Fig. 4, we can alter the in- 
teraction among an array of closely-packed resonators. 
In the usual case without the static magnetic field [Fig. 
4(a)], the electric field can transfer through the array 
due to the strong coupling between cylinders. It is be- 
cause each cylinder supports rotating states in opposite 
directions. When a clockwise rotational state is excited 
in one cylinder, the adjacent cylinder can attain the en- 
ergy by its state in anti-clockwise rotation. The process 
keeps on at other cylinders and a guiding of wave is pos- 
sible [shown as standing waves along the chain in Fig. 
4(a)]. Such a guiding effect of coupled waves has been 
discussed in details before 21 1. However, when we intro- 



duce the static magnetic field, only one rotational state 
is allowed in each cylinder at one frequency. When a 
clockwise rotational state is excited in one cylinder, the 
adjacent cylinder cannot attain the energy by its state in 
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First-order perturbation for Zeeman splitting 

Consider the Maxwell's equations: 
VxE = z'ft%fi-H 
VxB. = -icoe £-E 

For the sake of simplicity, we assume non-dispersive medium and construct the following eigenvalue 
problem: 



= co\ij/) , 
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ie £ Vx 




and iff) = 



It should be noted that the eigenvalue problem is subjected to certain boundary condition and the domain 
of interest. To avoid the problem of normalization for a state extending to infinity, we consider a finite 
circular domain with a radius several times larger than the radius of the ferrite cylinder and scattering 
boundary condition at the domain boundary. 



, we can reduce the 6x6 matrix equation 



to a 3x3 matrix equation. In cylindrical coordinates, the eigenvalue problem becomes 
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The eigenvalue problem can be written as 
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is the corresponding frequency operator for materials with isotropic permeability ju(r) and permittivity 
£ (r) , and 
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where is the unit step function. Note that for materials with isotropic permeability ju(r) and 
permittivity e(r) , the corresponding eigenvalue problem is 

^\yr) = ca\yf). 
Here, we have 
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For a system of cylindrical symmetry, the eigenstates can be written as | y/ n m 



H r (r) 
HAr) 



For small fl k <SC JU r , we can consider AQ as a first order perturbation to Q and divide the Zeeman 
effect into two steps. First, it shifts the frequency by replacing an isotropic permeability ju r to jLl' . 
Second, it splits the two degenerate eigenstates | y/ n m ^ and | \j/ n _ 
The splitting in first order perturbation theory is estimated as 



jjui l 6(a-r)\ jU oJ u(r)H;—^ + jU oJ u(r)H; 



(ys \An\Vi 



dr 



r dtp 



dA 



J(£ £(r)|E| 2 +// //(r)|H| 2 )^ 



2# r a ( .. ../rr* 
_ A) 



|(f £(r)|E| 2 +// //(r)|H| 2 )^ 



g r a f £ « ag, , ag 

2^- Ct) *°{ z dr dr 
A, |(£- £-(r)|E| 2 +// //(r)|H| 2 )^ 



2k 



m 



\E z (r = a)\ 2 



Ao \f y £ E(r)\E\ 2 +jU oJ u(r)\u\ 2 )dA 



In 



m£ ^\EXr = a)f 
2co n zK 



W jI( £o£(r )|E| 2 + MoM (r)\H\ 2 )dA 
Therefore, the frequency shift can be written as 
Ao) = 2kc 2 ■ ^ =a — 
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and C0 n is the corresponding eigen frequency of the unperturbed state. 
Here, we have used the usual definition of the angular momentum density: 

j = -^rx[E*x(VxE)]. 
2(Oi 

In our case, we have 
VxE = r^-f^ 



r d(p dr 

1 dp r)F 
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We can estimate the frequency shifts by evaluating U . As we mentioned before, there can be a problem 
of normalization when we consider an infinite spatial domain. To avoid this, we consider an approximate 
problem of finite cylindrical spatial domain of a radius several times larger than the ferrite cylinder. In 
this case, we can approximate U for this finite domain as 



= 2nE E d \ a Q E z 2 rdr 
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by using the fact that energy is mostly concentrated in the volume of the ferrite cylinder. 



By substituting E z (r) = E J m (k' n r) for an arbitrary constant E , we get 
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